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Abstract Young pulsars produce relativistic winds which 
interact with matter ejected during the supernova explo- 
sion and the surrounding interstellar gas. Particles are 
accelerated to very high energies somewhere in the pul- 
sar winds or at the shocks produced in collisions of the 
winds with the surrounding medium. As a result of in- 
teractions of relativistic leptons with the magnetic field 
and low energy radiation (of synchrotron origin, ther- 
mal, or microwave background), the non-thermal radia- 
tion is produced with the lowest possible energies up to 
~100 TeV. The high energy (TeV) 7-ray emission has 
been originally observed from the Crab Nebula and re- 
cently from several other objects. Recent observations 
by the HESS Cherenkov telescopes allow to study for the 
first time morphology of the sources of high energy emis- 
sion, showing unexpected spectral features. They might 
be also interpreted as due to acceleration of hadrons. 
However, theory of particle acceleration in the PWNc 
and models for production of radiation are still at their 
early stage of development since it becomes clear that 
realistic modeling of these objects should include their 
time evolution and three-dimensional geometry. In this 
paper we concentrate on the attempts to create a model 
for the high energy processes inside the PWNe which 
includes existence not only relativistic leptons but also 
hadrons inside the nebula. Such model should also take 
into account evolution of the nebula in time. Possible 
high energy expectations based on such a model are dis- 
cussed in the context of new observations. 
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1 Introduction 

The non-thermal nebulae around young energetic pul- 
sars (Pulsar Wind Nebulae - PWNe) have been sus- 
pected to accelerate leptons to sufficiently large ener- 
gies allowing production of TeV 7-rays (e.g. Gould 1965, 
Rieke & Weekes 1969, Grindlay & Hoffman 1971, Stepa- 
nian 1980). In fact, the nebula around the Crab pulsar 
was originally reported by the Whipple group as a first 
TeV 7-ray source (Weekes et al. 1989). Later, detections 
of TeV 7-rays from nebulae around other pulsars have 
also been claimed (e.g. around PSR 1706-44 - Kifune et 
al. 1995, Chadwick et al. 1998, Vela pulsar - Yoshikoshi et 
al. 1997, or PSR 1509-58 (MSH 15-52) - Sako et al. 2000), 
but some of the reported emission features have not been 
confirmed by recent more sensitive observations (see e.g. 
the results of HESS Collab. concerning PSR 1706-44 - 
Aharonian et al. 2005a, or Vela pulsar - Aharonian et 
al. 2006a). Recent advances in the study of the pulsar 
wind nebulae obtained mainly by observations in the 
lower energy range (from radio to X-ray emission) are 
reviewed by Gaensler & Slane (2006). 

Early detections of the TeV 7-ray emission from the 
PWNe have been usually interpreted in the so called 
synchrotron self-Compton model (SSC model) accord- 
ing to which relativistic leptons inside the nebula pro- 
duce soft synchrotron photons. These photons are next 
up-scattered by the same leptons to the 7-ray energies 
as a result of inverse Compton process (IC). Also soft 
photons of other nature, microwave background radia- 
tion (MBR), infrared or optical background present in- 
side the PNW, are up-scattered by these leptons to 7- 
ray energies. The semi-phenomenological model of this 
type has been elaborated by de Jager & Harding (1992) 
and successfully applied to the best studied TeV 7-ray 
source around the Crab pulsar. These authors apply this 
model assuming the distribution of the magnetic field 
inside the pulsar wind nebula obtained by Kennel & 
Coroniti (1984). Using the known spatial distribution of 
the synchrotron emission inside the nebula they get the 
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spectra and distribution of relativistic leptons inside the 
nebula. Next, from the obtained distribution of leptons 
and soft radiation inside the nebula, they calculate the 
TeV 7-ray emission (spectra and spatial distribution). 
Based on such procedure, de Jager & Harding predict 
that the higher energy TeV 7-ray emission should origi- 
nate closer to the pulsar. The model has been up-dated 
by Atoyan & Aharonian (1996) by including other pos- 
sible soft radiation targets inside the nebula and more 
recently confronted with the observations of the Crab 
Nebula by Hillas et al. (1998). For the application of SSC 
model to nebulae around other pulsars see e.g. du Plessis 
et al. (1995) - PSR 1509-58, or de Jager et al. (1996) - 
Vela pulsar. 

The classical model of de Jager & Harding (1992) 
applies the multiwavelength observations of the Crab 
Nebula predicting the TeV 7-ray spectra at the present 
time. It does not consider the evolution of the nebula 
in time. The first simple time dependent leptonic model 
for the PWNe has been discussed by Aharonian, Atoyan 
& Kifune (1997) and applied to the nebula around PSR 
1706-44. The authors analyze the evolution of the equi- 
librium energy spectrum of leptons injected at a constant 
rate, with a fixed spectrum (independent of time) into 
the region with a constant magnetic field. They come 
to an interesting conclusion that PWNe with relatively 
low magnetic field should become strong sources of 7-ray 
emission. Therefore, the efficiency of energy conversion 
from the rotating pulsar to 7-rays should increase with 
the age of the nebula. 

Observation of the TeV 7-ray emission up to ^80 
TeV by the HEGRA Collab. (Aharonian et al. 2004) 
indicates the existence of particles with energies up to 
~10 15 eV inside the Crab Nebula. Acceleration of lep- 
tons to such energies requires special conditions in the 
acceleration region (low synchrotron energy losses, very 
efficient acceleration mechanism). Note however, that re- 
cently the HESS Collab. (Aharonian et al. 2006b) re- 
ported the 7-ray spectrum from the Crab Nebula show- 
ing an exponential cut-off at energy ~14 TeV. If real 
such discrepancy might suggest the presence of an ad- 
ditional component in the Crab Nebula spectrum above 
several TeV, which could be interpreted as the contribu- 
tion from relativistic hadrons inside the nebula (see e.g. 
Bednarek & Bartosik 2003, BB03). In fact, the existence 
of relativistic hadrons inside the PWNe has been con- 
sidered in the past by e.g., Cheng et al. (1990), Atoyan 
& Aharonian (1996), and Bednarek & Protheroe (1997). 
These hadrons interact with the matter of the super- 
nova remnant and surrounding medium producing pions 
which decay to 7-rays and neutrinos. Possible contribu- 
tion of 7-rays from hadronic processes to the observed 
flux from the Crab Nebula has been discussed in the 
above mentioned papers. For example, Atoyan & Aha- 
ronian (1996) consider interaction of relativistic leptons 
and hadrons with the matter inside the Crab nebula and 
conclude that bremsstrahlung and tt° processes might 



give some interesting contribution to the observed spec- 
trum provided that the effective density of matter inside 
the nebula is by one order of magnitude larger than the 
observed average density of matter. In this paper we dis- 
cuss a more complete hadronic-leptonic model which has 
been recently proposed independently by Bednarek & 
Bartosik (2003) and Amato, Guetta & Blasi (2003). Let 
us first review some new observational results in which 
context this model will be considered. 



2 New observations of the PWNe 

Since the aim of this paper is to review the expected 
multi-messenger high energy signatures for acceleration 
of particles in the PWNe, here we only mention the most 
interesting (in our opinion) new observational results ob- 
tained mainly by the HESS Collaboration in the TeV 
7-ray energy range. For a more complete review of the 
HESS results see e.g., de Jager (2006), Gallant (2006), 
or Aharonian (2006). 

The great advantage of the HESS Cherenkov tele- 
scope system is the ability of morphological studies of 
TeV 7-ray emission. It has appeared that TeV 7-ray 
emission from the middle-aged PWNe is extended and 
comes from the region up to several parsecs, e.g. MSH 
15-52 (Aharonian et al. 2005b). Moreover, the emission 
is often asymmetric which is probably the result of inter- 
action of the supernova with the surrounding medium, 
e.g. PWNe associated with PSR B1823-13 (Aharonian et 
al. 2005c), Kookaburra complex (Aharonian et al. 2006c) 
or around the Vela pulsar (Aharonian et al. 2006a). Such 
asymmetries are difficult to explain by the movement of 
the pulsars through the nebulae (e.g. the observed off- 
set of the Vela pulsar by ^0.5 degree would require its 
velocity of ^250 km s _1 , while the measured value is 
only 65 km s _1 , Caraveo et al. 2001). Another interest- 
ing result of the HESS Collab. is the discovery of 7-ray 
emission between 550 GeV up to 65 TeV from the region 
of the Vela pulsar with unexpected very flat spectrum 
(index ^1.45) and the cut-off at ^14 TeV (Aharonian 
et al. 2006a). This extended TeV emission offset from 
the location of the Vela pulsar has been also confirmed 
the CANGAROO Collab. (Enomoto et al. 2006). In fact, 
all other PWNe have the spectral indices in the range 
~2.2-2.4, i.e. more similar to that observed for the Crab 
Nebula. The TeV 7-ray spectral features observed from 
the Vela nebula seem to be more consistent with their 
hadronic origin (e.g. Bednarek & Bartosik 2003, Horns 
et al. 2006). 

Since in the last few years no essential development 
in the modeling of the high energy emission from the 
PWNe in the framework of the pure leptonic models has 
been obtained, we concentrate here on the recently con- 
sidered hadronic-leptonic (hybrid) models which try to 
include not only radiation processes due to the injected 
leptons but also processes due to injection of relativistic 
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hadrons into the nebula. These recent works were stimu- 
lated mainly by the studies of the mechanism of possible 
acceleration of leptons as a result of their interaction 
with the Alfven waves generated by coherently gyrating 
heavy nuclei in the PWNe (for review see Arons 1998). 
Such hybrid lcptonic-hadronic models also predict addi- 
tional high energy messengers (neutrinos, neutrons or 
charged particles contributing to the observed cosmic 
rays at the Earth), whose observation will certainly put 
new constraints on the high energy processes occurring 
in the PWNe. 



3 Hadronic-leptonic (hybrid) model 

Possible acceleration of heavy nuclei, extracted from the 
surface of the neutron star, by the electric field induced 
in the outer gaps of the inner pulsar magnetosphere was 
suggested already by Cheng et al. (1986). However, these 
nuclei should partially photo-disintegrate when passing 
through the non-thermal radiation of the outer gaps. 
Neutrons from their disintegration are injected into ex- 
panding nebula surrounding the energetic pulsar. The 
decay products of neutrons, i.e. protons, partially cap- 
tured inside the nebula, interact with the matter and 
contribute to the 7-ray emission. Such model was sug- 
gested for the 7-ray and neutrino production from the 
Crab Nebula (Bednarek & Protheroe 1997). 

The nuclei can be additionally accelerated by the pul- 
sar wind carrying significant amount of energy lost by 
the pulsar. They can energetically dominate the rela- 
tivistic pulsar wind as proposed by Hoshino et al. (1992) 
and Gallant & Arons (1994). Based on the observations 
of the wisps in the Crab Nebula, these authors conclude 
that the nuclei are accelerated somewhere inside the pul- 
sar wind zone (without specifying the details of the un- 
known acceleration process). After passing the pulsar 
wind shock, the nuclei generate Alfven waves which en- 
ergy is next absorbed by leptons. As a result, monoen- 
ergetic nuclei and leptons with a power law spectrum 
are injected into the nebula. The radiation model for 
the PWNe based on the above mentioned acceleration 
mechanism of leptons has been independently developed 
by Amato, Guetta & Blasi (2003) and Bednarek & Bar- 
tosik (2003). These models try to include (for the first 
time self-consistently) relevant leptonic and hadronic ra- 
diation processes inside the nebula taking into account 
also the time evolution of the nebula. We describe the 
main features of such a model and obtain results based on 
the approach presented in Bednarek & Bartosik. Note, 
that although both models differ in some details, they 
give not very different results in the case of the Crab 
Nebula the concerning contribution of hadronic processes 
to the highest energy part of the 7-ray spectrum and 
possible detection of the pulsar wind nebulae by the fu- 
ture neutrino telescopes. Such hybrid hadronic-leptonic 
model has been also considered as possible explanation 



of the 7-ray emission from the recently discovered 7- 
ray nebula associated with the Vela pulsar (Horns et 
al. 2006). 

3.1 A simple model for expansion of the nebula 

The time evolution of a supernova remnant under the 
influence of an energetic pulsar is described according to 
the general picture proposed by Ostriker & Gunn (1971) 
and Rees & Gunn (1974). Let us denote the initial expan- 
sion velocity of the bulk matter in supernova envelope by 
Vo.sn and its initial mass by Mo,sn- The expansion ve- 
locity can increase due to additional supply of energy to 
the nebula by the pulsar. It can also decrease due to ac- 
cumulation of the surrounding matter. These processes 
are taken into account in order to determine the radius 
of the nebula at a specific time, t, by using the energy 
conservation, 

^^ = ^ + fM^ (i) 
2 2 j 

where 

L pul (t) = B 2 s Rf^/6c 3 « 3 x 10 43 B2 2 P m 4 erg s -i ; (2 ) 

is the pulsar energy loss, Q = 2n/P, and P = 10 _3 P ms 
s changes with time according to P^ s (t) — Po yms + 2 x 
10~ 9 tBi2, Po,ms is the initial period of the pulsar, and 
B = 10 12 i?i2 G is the surface magnetic field of the pulsar. 
The nebula increases its mass according to 

M SN (*) = M , SN + -TTp SUI R^ eh {t), (3) 

where p sur is the density of surrounding medium and 
i?Ncb is the radius of the expanding envelope at the time 

t, 

i?Ncb- / V SN (t')dt'. (4) 

Jo 

The average density of matter inside the nebula is 

PNcb = 3M SN (t)/47ri& eb (t). (5) 

The above set of equations is solved numerically in order 
to obtain the basic parameters of the expanding PWN. 

Also the location of the inner shock inside the PWN 
can be estimated by balancing the pressure of the ex- 
panding nebula with the pressure of the pulsar wind 
(Rees & Gunn 1974), 

Vg) ti<rL^{t')dt' 

where a is the ratio of the magnetic energy flux to the 
total energy flux lost by the pulsar at the location of 
the pulsar wind shock at radius R s h- Then, the magnetic 
field strength at the shock region can be estimated from 
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where i? pu i and £> pu i are the radius and the surface mag- 
netic field of the pulsar, and R\ c is the light cylinder 
radius. 

The evolution of basic parameters characterizing the 
pulsar wind nebula, for some initial parameters of the 
supernova explosion and the pulsar, are shown on Fig. 1 
in BB03. 

3.2 Model for acceleration of leptons and hadrons 

As already noted above we are interested in the scenario 
in which rotating magnetospheres of neutron stars can 
accelerate not only leptons but also heavy nuclei, ex- 
tracted from positively charged polar cap regions. In fact, 
different aspects of the high energy phenomena around 
pulsars, such as the change in the drift direction of the 
radio sub-pulses (Gil et al. 2003), the existence of mor- 
phological features inside the Crab Nebula (so called ra- 
dio wisps), and the appearance of extremely energetic 
leptons inside it (Gallant & Arons 1994), can be natu- 
rally explained by the presence of heavy nuclei. Arons 
and collaborators (e.g. see Arons 1998) postulate that 
the Lorentz factors of iron nuclei accelerated somewhere 
in the inner magnetosphere and/or the Crab pulsar wind 
zone should be, 

7Fe « riZe$ opcn /m Fc c 2 w 8 x 10 9 77-Bi 2 P~ s 2 , (8) 

where TOFe and Ze are the mass and charge of t he iron 
nuclei, c is the velocity of light, and $ op cn = \/L pu \/c 
is the total electric potential drop across the open mag- 
netosphere, and r] is the acceleration factor determin- 
ing the Lorentz factor of nuclei in respect to the max- 
imum one allowed by the pulsar electrodynamics. Fol- 
lowing Arons and collaborators, the authors assumed: 
(1) r/ is not very far from unity, the value r\ = 0.5 is 
adopted; (2) iron nuclei take most of the spin down power 
of the pulsar, Lp c = xL mt , where \ = 0.95. Unfortu- 
nately, this values are not predicted at present by any 
model of the ion acceleration in the pulsar wind and can 
only be constrained by the high energy observations of 
the PWNe. The iron nuclei are extracted from the neu- 
tron star surface and accelerated during the pulsar radio 
phase when the efficient leptonic cascades heat the polar 
cup region. They are farther accelerated in the pulsar 
wind zone due to the linear increase of the Lorentz fac- 
tor of the outflowing plasma (as more recently discussed 
by Contopoulos & Kazanas 2002) or due to reconnec- 
tion of oppositely directed magnetic fields in the wind 
(Michel 1982, Coroniti 1990, Lyubarsky & Kirk 2001). 
The quasi-monocncrgctic heavy nuclei, after crossing the 
pulsar wind shock, generate Alfven waves in the down- 
stream region, which energy is resonantly transfcrcd to 
leptons present in the wind (Hoshino et al. 1992). As 
a result, leptons are accelerated to high energies with 
a power law spectrum with the spectral index b\ s=s 2 
between E\ = 7p e m e c 2 and E 2 w ■yFeAm p c 2 /Z (see Gal- 
lant & Arons 1994), where m c and m p are the electron 



and proton mass, respectively. The spectrum is normal- 
ized in such a way as to get the conversion efficiency 
of energy from the iron nuclei to the leptons equal to 
£. Note that the radiation from leptons depends on the 
product of the energy conversion from the pulsar wind to 
nuclei, x, and the acceleration efficiency of positrons by 
these ions, £. Therefore, decreasing the first coefficient 
and increasing the second, one obtains the same level of 
radiation from positrons but a lower level of gamma-ray 
flux from hadronic interactions of ions with the matter 
inside the nebula. Since the dependence of % • £ on time 
is not predicted by any theoretical model in the original 
paper by Bednarck & Bartosik (BB03), this value was 
kept constant during the evolution of the nebula. Rcla- 
tivistic particles accelerated by the mechanism discussed 
above are captured inside the pulsar wind nebula losing 
energy on different processes. However, nuclei injected 
into the nebula at a specific time i; n j, escape from it at 
the time t osc , if their diffusion distance in the turbulent 
magnetic field of the nebula, -Rdiff , is equal to the dimen- 
sion of the nebula, i?Nob, at the time t csc . The diffusion 
distance in the magnetic field of the nebula is obtained 
by integration (Bcdnarek & Protheroe 2002), 

RdiS = / V w*'' (9) 

where the diffusion coefficient is taken to be D = R^c/3, 
and i?L is the Larmor radius of nuclei depending on the 
distance from the center of the nebula. i?diff and t csc 
are calculated numerically. The adiabatic losses of nu- 
clei during their propagation inside the nebula are also 
included. Due to these losses and interactions of nuclei 
with the matter the energy of a nucleus at a specific time 
is 

B (*) = ^)|r^> ( 10 ) 

where -E(tj n j) is its energy at the time fj n j, K is the inelas- 
ticity coefficient for collisions of nuclei with matter, and 
r is the optical depth calculated from the known density 
of matter inside the nebula. The details of these calcu- 
lations are given in sections 4.2 and 4.3 in Bcdnarek & 
Protheroe (2002). All the above mentioned effects have 
been taken into account when considering leptons and 
hadrons inside the nebula. 



4 High energy radiation from the PWNe in the 
hybrid model 

Having defined the model of the nebula expansion (which 
gives the basic parameters of the nebula), and also the 
injection spectra of particles as a function of time after 
the supernova explosion (the age of the nebula) , the equi- 
librium spectra of leptons and hadrons inside the nebula 
can be obtained taking into account different energy loss 
processes. Leptons injected into the medium of the ex- 
panding supernova remnant suffer energy losses mainly 
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on radiation processes, bremsstrahlung, synchrotron, and 
the inverse Compton, and due to the expansion of the 
nebula. The rate of their energy losses can be described 

by 

dE 

- -77- = (ai + a 2 )E + (ft + p 2 )E 2 GeV s"\ (11) 
at 

where a.\ and a 2 describe the bremsstrahlung and adi- 
abatic losses; [3\ and fa the synchrotron and ICS en- 
ergy losses. The coefficients, a\,a 2 , Pi, and [3 2 , depend 
on time in a complicated way due to the changing con- 
ditions in the expanding nebula (magnetic field, density 
of matter and radiation). Therefore, the above formula 
can not be integrated analytically at an arbitrary time 
after supernova explosion. In order to determine energies 
of leptons, E, inside the nebula at a specific time t ohs , 
which have been injected with energies E a at an earlier 
time t, we use the numerical approach. However, know- 
ing the parameters of the nebula determined at a given 
time t, the evolution of the equilibrium spectrum of lep- 
tons during the time step At is determined analytically. 
Next, the conditions inside the nebula are changed to 
values which are obtained from the expansion model of 
the nebula at time t + At. The equilibrium spectrum of 
leptons at, i Q b s , is then obtained by summing over the 
spectra injected at specific time and over all time steps 
up to the present observed time t t, s , 



dN(t ohs ) 
dE 



dN 
dE~dt 



■dt. 



(12) 



where dN/dE a dt is the injection spectrum of leptons at 
time t, t' = t bs — t, and the jacobian J(t') = E a /E(t') 
describes the change of energy of lepton during the pe- 
riod At. It is calculated analytically by solving Eq. 11. 
The example equilibrium spectra of leptons inside the 
nebula at the specific time after explosion of supernova 
are given in Fig. 3 in BB03. 

In order to calculate the equilibrium spectra of dif- 
ferent types of nuclei inside the nebula at its given age, 
we take into account the energy losses on the interaction 
of nuclei with the matter of the nebula, the adiabatic en- 
ergy losses, and the escape of nuclei from the nebula due 
to diffusion. The example equilibrium spectra of nuclei 
inside the nebula at specific age are shown in Fig. 2 in 
Bednarek & Bartosik (2003). 



4.1 Gamma-rays 

Having obtained the equilibrium spectra of relativistic 
leptons and hadrons as a function of time after super- 
nova explosion we can calculate the energy spectra of 
photons produced inside the nebula by these particles 
in different radiation processes. In general, such spectra 
are characterized by three main components: two due to 
leptons (synchrotron, IC) and one due to hadrons (n° 
decay). Depending on the parameters of the pulsar, the 
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Fig. 1 The 7-ray spectrum from Vela nebula modeled with 
hybrid hadronic-leptonic model (Bednarek & Bartosik 2003) . 
7-rays from 7r° decay created in hadronic collisions (dot- 
dashed curve) and from ICS of leptons (full curves). The 
acceleration mechanism of leptons operates through all the 
lifetime of the pulsar (the thinnest full curve), or only during 
the first 5000 yrs, 3000 yrs, and 2000 yrs (the thickest full 
curve) . The acceleration of hadrons is continuous throughout 
the age of the pulsar ~11000 yrs. The dashed and dotted 
curves show the synchrotron emission from the extended and 
compact nebulae. 



supernova, the surrounding medium and on the age of 
the nebula, different components may dominate. For the 
PWNe in a very early stage of development (less than 
^100 yrs), 7-rays from the interaction of hadrons with 
matter may dominate. For young nebulae (at the age of 
~10 3 yrs, the Crab type), the power of the synchrotron 
spectrum dominates over that of the IC 7-ray spectrum. 
In the case of the middle age nebulae (at the age ^10 4 
yrs, the Vela type), the power in synchrotron bump may 
become comparable to the power in the TeV 7-ray bump 
(produced mainly due to the IC scattering of the MBR). 

For reasonable initial parameters, the hybrid model 
predicts also the appearance of the small bump in the 
multiwavclcngth spectrum of the Crab nebula at the 
highest energies, see Bednarek & Bartosik (2003, 2005) 
or Amato, Guetta & Blasi (2003). In fact, recent mea- 
surements of the cutoff in the Crab Nebula spectrum 
at <~14 TcV (Aharonian et al. 2006a), combined with 
the previous measurements of the spectrum up to ^80 
TeV (well described by a simple power law, Aharonian 
et al. 2004), might be interpreted as the evidence of such 
an additional component at ^100 TeV. 

Based on the hybrid model, Bednarek & Bartosik 
(BB03) calculated expected emission from the nebula 
around the Vela pulsar predicting significant contribu- 
tion to the 7-ray flux from the decay of n°. However, in 
the original calculations it has been assumed that rela- 
tivistic leptons are injected into the nebula as a result of 
resonant interactions with nuclei during the whole life- 
time of the pulsar. Recent results obtained by the HESS 
Collab. (Aharonian et al. 2006a) show that it may not 
be the case. As we noted above, the observed TeV 7-ray 
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spectrum from the Vela Nebula is flatter than E~ 2 , and 
very nicely resembles the pure hadronic component (see 
Fig. 7 in BB03). Fig. 1 shows the 7-ray spectra obtained 
in terms of the hybrid model which are re-calculated ap- 
plying the modified assumption about the acceleration 
process of leptons by hadrons inside the nebula. In the 
present calculations, leptons are accelerated only during 
the first 2000, 3000, and 5000 yrs after the pulsar birth. 
Leptons accelerated in the past lose efficiently their en- 
ergy on radiative and adiabatic processes. Therefore, the 
equilibrium spectrum of leptons inside the nebula shifts 
to lower energies with increasing age of the nebula due 
to the lack of freshly injected leptons at the present time. 
As a result, the IC spectrum obtained in such a modi- 
fied model cuts-off at lower energies and the flat 7-ray 
emission above ~1 TcV due to the decay of ir° can be 
clearly observed. Such model predicts the appearance of 
a steep component in the TeV 7-ray spectrum below ~1 
TeV which might be used to estimate the duration of the 
acceleration process of leptons. 

Due to the proximity of the pulsar and diffusion of 
leptons from the acceleration site close to the pulsar wind 
shock, it has been predicted that the TeV 7-ray emission 
from the Vela PWN should be extended with charac- 
teristic dimension of the order of ~5 pc (Bednarek & 
Bartosik 2005). This is consistent with the HESS obser- 
vations. Note, that the Vela nebula is significantly dis- 
placed from the location of the Vela pulsar which might 
be the result of unequal pressure on the nebula due to 
its interaction with the inhomogeneous medium as ear- 
lier suggested by Blondin, Chavalier & Frierson (2001). 

Note, however that the morphological studies of the 
TeV 7-ray emission features of the pulsar wind nebula 
G18.0-0.7, close to the pulsar PSR J1825-137, show clear 
softening of the spectrum with increasing distance from 
the pulsar (Aharonian et al. 2006d). Such feature is bet- 
ter explained by the leptonic model in which leptons lose 
radiatively and adiabatically energy when diffusing out- 
side the pulsar. 



4.2 Neutrinos and Neutrons 

Provided that the hypothesis on hadron acceleration in 
some PWNe is correct and assuming that significant 
amount of observed 7-ray flux is due to hadrons, it is 
possible to estimate the expected neutrino flux and the 
event rates in the presently planned and constructed 1 
km 2 neutrino detectors (IceCube, KM 3 NET). Such sim- 
ple estimate, based on the observed TeV 7-ray flux, has 
been performed by Amato & Guetta (2003) for a few 
PWNe. For example these authors argue that ~12 neu- 
trino events per km 3 per year should be detected from 
the Crab Nebula. More detailed calculations according 
to the hybrid model of Amato, Guetta & Blasi (2003) 
predicts ~5-13 neutrino events per yr per km 3 from the 
Crab Nebula, provided that hadron Lorentz factors reach 



~10 4 -10 7 . The event rates estimated for the Crab and 
Vela pulsars, based on a similar model by Bednarek (2003), 
are on the level of ~1 neutrino event within 1 km 2 de- 
tector during 1 year, i.e at the level of the atmospheric 
neutrino background (ANB). However, when these nebu- 
lae were younger, then the event rates were significantly 
larger, decreasing with the age of the nebula, i n cb, pro- 
portionally to ~ t~ eh (see Fig. I in Bednarck 2003). 
Based on the observed fluxes of the TeV 7-rays from 
the PWNe by the HESS Collab., Kappes et al. (2006) 
predicts on average 1.4 — 4.1 neutrino events per year 
above >5 TeV from the direction of the Vela Nebula in 
the KM 3 NET detector, whereas the estimated ANB is 
of the level of ~1 neutrino event. 

The hybrid model for 7-ray production in the PWNe 
postulates also that from relatively close sources, the 7- 
rays and neutrinos might be accompanied by neutrons. 
In fact, 7-ray emission up to ~65 TeV from the direction 
of the Vela Nebula, requires the acceleration of hadrons 
to the Lorentz factors of the order of ~10 6 . In the past, 
when the Vela pulsar was at the present age of the Crab 
pulsar, nuclei might be accelerated to the Lorentz factors 
of the order of magnitude larger, i.e. ^10 7 . The mean 
free path for neutrons extracted from nuclei with such 
Lorentz factors (or produced in hadronic interactions of 
protons) is already ~100 pc. So then, significant number 
of neutrons from the Vela pulsar might also reach the 
Earth about 10 4 yrs ago, from the distance of <~ 300 pc 
(Caraveo et al. 2001). 

4.3 Contribution to Cosmic rays? 

Accepting the hybrid model and keeping in mind that 
PWNe are at present the best established TcV 7-ray 
sources, we should expect that particles escaping from 
the PWNe might contribute to the cosmic rays (CR) 
in the Galaxy. In fact, pulsars have been suspected since 
their discovery as main sources of cosmic rays. Ostriker & 
Gunn (1969) and Karakula, Osborne & Wdowczyk (1974) 
postulated acceleration of CRs above the knee region in 
the pulsar wind zones by the large amplitude electro- 
magnetic waves generated by rotating neutron stars. The 
contribution of particles accelerated by pulsars to the 
observed cosmic ray spectrum has was later discussed 
by Cheng & Chi (1996), Bednarek & Protheroe (2002), 
Giller & Lipski (2002). Cheng & Chi (1996) propose that 
nuclei injected from the pulsar magnetospheres during 
the pulsar glitches can contribute to the knee region in 
the CR spectrum. Ciller & Lipski (2002) derive the initial 
parameters of the pulsar population inside the Galaxy re- 
quired to explain the observed shape of the CR spectrum 
and its intensity up to the highest energies. Bednarck & 
Protheroe (2002) estimate the contribution of heavy nu- 
clei accelerated in the pulsar outer gaps (Cheng, Ho & 
Ruderman 1986) from the population of Galactic radio 
pulsars to the CRs above the knee region taking into ac- 
count the propagation and escape conditions, adiabatic 
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and collisional energy losses of the nuclei during their 
propagation in the pulsar wind nebulae. In fact, a new 
component in the CR spectrum at energies above <~ 10 15 
eV, supplied e.g. by pulsars, is required by the measure- 
ments of the mass composition, which suggests an in- 
crease of the average mass above the knee (e.g. Glas- 
macher et al. 1999, Ave et al. 2003). 

Based on the hybrid model, Bednarek & Bartosik 
(2004) calculated the spectra of nuclei which escape from 
the PWNc around the pulsars with assumed initial pa- 
rameters. Applying a few different models for the initial 
parameters of the pulsars formed inside the Galaxy, con- 
cerning the distribution of the surface magnetic fields of 
the new born neutron stars (derived by Narayan 1987) 
and the initial pulsar periods (derived by Lorimer et 
al. 1993), the contribution of the pulsar population to 
the CRs in the Galaxy has been estimated. It is con- 
cluded that model B for the pulsar population presented 
in Lorimer at al. (1993) is able to describe satisfactory 
the shape of the CR spectrum and the features of the 
mass composition in the energy range between the knee 
and the ankle (see Figs. 3b and 4 in Bednarek & Bar- 
tosik 2004). The fine features in the CR spectrum (the 
knee region) might also appear due to the presence of the 
nearby pulsars (e.g. Erlykin & Wolfendale 2004). How- 
ever, such pulsars would create an anisotropy above the 
present observational limit (Bhadra 2006). 

Also extremely high energy CRs, i.e. above the ankle 
at <~ 10 18 eV, might be accelerated in the wind regions 
of pulsars with super strong surface magnetic fields, so 
called magnetars. Blasi, Epstein & Olinto (2000), pro- 
pose that the highest energy CRs are accelerated locally 
in our Galaxy by the pulsars whose initial periods are 
shorter than ~10 ms and whose surface magnetic fields 
are 10 12 — 10 14 G. If such pulsars inject heavy nuclei with 
energies ~ 10 20 eV, then their Larmor radii are smaller 
than the dimensions of the galactic halo, allowing their 
efficient trapping inside the Galaxy. The authors argue 
that such nuclei can be accelerated to extremely high 
energies after the moment when supernova envelope be- 
comes transparent. However it is not clear if they are not 
captured inside the supernova envelope by the magnetic 
field of the PWN? For example, from the modeling of the 
PWNe (BB03), it is possible to estimate the magnetic 
field inside the nebula, ~10-100 G, and its dimension, 
<~ 10 15 — 10 16 cm, at ~1 year after explosion. For such 
parameters, the Larmor radius for iron nuclei are lower 
than the dimensions of the nebula allowing an efficient 
trapping. Then, the adiabatic energy losses of nuclei due 
to fast expansion of the nebula become important. In a 
similar model, Arons (2003) suggests that pulsars with 
extreme parameters produced in the whole Universe are 
responsible for the highest energy CRs. The winds of 
such pulsars disrupt the supernova envelopes allowing in 
this way acceleration of protons to extremely high ener- 
gies (of the order of 10 21 ~ 22 eV), and their escape from 
the nebula. The final shape of the CR spectrum at the 



highest energies in the universe depends on the rate of 
gravitational energy losses by the pulsar during its early 
stage after formation. 

Note that the PWNe associated with such very young 
pulsars in the Galaxy should also produce significant 
fluxes of neutrinos which might be detectable within a 
few years after their birth, see e.g. Berezinsky & Prilut- 
sky (1978), Protheroe, Bednarek & Luo (1998), Beall & 
Bednarek (2002). 



5 Interaction of PWNe with massive clouds 

As reported in Sect. 2, the TeV sources are often sub- 
stantially displaced from the positions of the associated 
middle aged pulsars. This is probably the effect of the 
interaction of a parent supernova shock wave with the 
surrounding inhomogeneous medium. This expectation 
seems to be very natural since young supernovae are usu- 
ally close to the star forming regions containing massive 
clouds. Particles escaping from the PWNe can be par- 
tially captured by these high density clouds producing 
7-rays in collisions with the cloud matter. In fact, the 
TeV 7-ray emission from one of the most massive com- 
pact stellar association, Cyg OB2, has been detected by 
HEGRA Collab. (Aharonian et al. 2002). Observation of 
other 14 young open clusters does not give positive de- 
tections due to the insufficient sensitivity (Aharonian et 
al. 2006e) . Recently extended TeV 7-ray emission, coin- 
cident with the location of giant molecular clouds, has 
been also observed from the direction of the Galactic 
Center (Aharonian et al. 2006f). Such extended 7-ray 
sources can be also caused by the appearance of nearby 
energetic pulsars. If such pulsars have been born typi- 
cally ~10 4 yrs ago with initial periods of the order of 
a few millisecond and the surface magnetic fields of the 
order of ~10 13 G, then at present they should have peri- 
ods close to ^200 ms. Such radio pulsars are difficult to 
discover in the high density regions. Hadrons accelerated 
during the lifetime of a pulsar can be partially captured 
in dense (10 3 - 10 5 cm~ 3 ), magnetized (10" 5 - 10~ 3 G) 
molecular clouds, producing 7-rays for a relatively long 
time. Note, that even hadrons with Lorentz factors up 
to ^10 9 might be captured in the molecular clouds with 
dimensions of several parsecs since their Larmor radii are 
small enough. Possible radiation effects from the interac- 
tions of the PWNe with high density regions have been 
considered in the context of the 7-ray emission from the 
Galactic Center (Bednarek 2002), and also proposed as a 
possible explanation of the unidentified TeV 7-ray source 
in the Cyg OB2 region (Bednarek 2003). 



6 Conclusion 

Pulsars produce a class of the best established TeV 7-ray 
sources in which particles are accelerated above ^100 
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TeV. Observations with the HESS telescopes allow for 
the first time to study morphology of these objects in- 
dicating that the interaction of parent supernova rem- 
nants with the surrounding medium strongly influences 
the evolution of the pulsar wind nebulae. These interac- 
tion effects have not been considered in the past leptonic 
models for the production of high energy radiation inside 
PWNe and strongly encourage the efforts to work on a 
much more advanced, time dependent, three dimensional 
model for the PWNe. Moreover, new observations of the 
PWNe at TeV 7-rays seem to suggest that acceleration of 
hadrons to high energies inside the PWNe should be also 
seriously considered. In fact, the unusual shape of the 7- 
ray spectrum from the nebula associated with the Vela 
pulsar (relatively narrow peak at several TeV) and possi- 
bly weak evidence of a discrepancy between observations 
of the Crab Nebula at the highest energies by HEGRA 
and HESS can be naturally explained by the existence of 
an additional component in the TeV 7-ray spectra from 
these two objects which is due to contribution from rel- 
ativistic hadrons. Since hadrons are weakly coupled to 
matter, they can diffuse from the PWNe, interact with 
the surrounding medium, and also contribute to the ob- 
served cosmic ray spectrum in the Galaxy. Therefore, 
PWNe should be considered as one of the serious candi- 
dates responsible for the bulk of the cosmic rays above 
the knee region, i.e. >10 15 eV in the Galaxy. Detection 
(or non-detection) of other neutral radiation from the 
PWNe, i.e. neutrinos and neutrons, should give the an- 
swer to this question. 

The aim of this paper was to stress the importance 
of such multi-messenger high energy observations of pro- 
cesses occurring inside the PWNe. 
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